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Abstract: We study the classical stability of a class of S-brane geometries having 
cosmological horizons. By considering the perturbations of the metric in these ge- 
D . ometries we establish that their horizons are unstable in the sense that an observer 



trying to cross the horizon experiences an infinite flux of radiation at the instant 
of crossing. The backreaction of this radiation is likely to convert the horizons into 
curvature singularities, similar to the instability of the internal Cauchy horizon of 
the Reissner-Nordstrom black hole. We also compute the particle production by 
the time-dependent fields in the future regions of these geometries, and find that 
the spectrum of produced particles is thermal, with temperature coinciding with the 
Hawking temperature computed by euclideanizing the metric in the static region. 
Possible implications of these results are discussed. 
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1. Introduction 

S-branes are spacelike surfaces in spacetime along which transitions between different 
vacua are speculated to take place in string theory [p]]. Their study is motivated by 
the desire to be able to extend the present theoretical toolbag to include techniques 
for understanding time- dependent problems in string theory. If the transitions de- 
scribed by S-branes really do arise, parts of their behavior should be describable in 
terms of the evolution of the low-energy fields of supergravity, perhaps also with a 
rolling tachyon or tachyons which describe the field-theory version of the transition 
between vacua. 

Several time-dependent supergravity solutions have been proposed as corre- 
sponding to S-branes [§, || |], [5], §]. Although many have an FRW-like singularity at 
early times, some do not and instead have an interesting global structure consisting 
of non-singular and asymptotically-flat past and future regions separated by static 
regions having time-like singularities. Some of these geometries have been studied in 
fair detail. In particular, Ref. || performs a general analysis of the charge, tension, 
entropy and Hawking temperature of many of these spaces. It was concluded there 
that the time-dependent regions provide an interesting interpretation that fits very 
well with the original S-brane proposal H . The static regions could be interpreted as 
the fields external to a pair of oppositely-charged, negative-tension branes. A similar 
interpretation was proposed earlier in terms of orientifold planes (7]. 

Our purpose here is to explore two important issues of stability for these ge- 
ometries, which were not completely addressed in Ref. 0. The first of these is the 
classical stability of the solutions. Our preliminary study showed that small per- 
turbations in Klein-Gordon scalar fields in these geometries do not grow with time, 
but do infinitely blue-shift for inertial observers passing through the horizons. This 
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signals a potential instability, because the backreaction of this energy density on 
the metric is ultimately likely to convert the horizons into curvature singularities 
RJ. Here we extend this observation somewhat by performing a similar analysis 



for the modes of the metric itself, adapting for this purpose the methods applied by 
Chandrasekhar and Hartle to the Reissner-Nordstrom black hole. We arrive at the 
same conclusion as before, thereby strengthening the arguments that the horizons of 
this geometry are unstable. This fact indicates that gravitational perturbation will 
generally introduce null-like singularities in the geometry. 

The presence of these singularities makes these solutions more similar in char- 
acter to the general S-brane solutions discussed in ||, for which similar curvature 
singularities are already present on null-like surfaces at the classical level. Only for 
special values of some parameters do the singular surfaces of || become horizons, 
and based on our calculation one might wonder if these horizons are stable to the 
formation of singularities once metric fluctuations are considered. 

The second question we address concerns the particle production which is pro- 
duced by the non-static metric at late times. In ref. M a Hawking temperature was 
associated with the static regions of the metric, by requiring in the usual way that 
there be no conical singularities in their Euclidean sections. This temperature was 
argued to be related to the particle spectrum which is seen by the (accelerating) 
static observers. Here we directly compute the particle production in the non-static 
region and show that it is also thermal in character, with the same temperature that 
was found previously for the static regions. 



2. Classical Instability 

In this section we present in detail the stability analysis for gravitational metric 
perturbations of the simplest SO-brane geometry of ref. 0. We will closely follow a 
procedure used by Chandrasekhar and Hartle [§, ||, to show the instability of the 
Reissner-Nordstrom (RN) metric. We show that the same conclusion also applies to 
the SO-brane geometries of interest here. 

We start by reviewing the preliminary analysis performed in ref. where the 
Klein-Gordon modes were analysed arriving to the preliminary conclusion of instabil- 
ity of the Cauchy horizon. We then perform a more rigorous analysis by considering 
in detail the metric perturbation modes, arriving to the same conclusion of instability. 

Consider then the simplest geometry that corresponds to an uncharged SO-brane 
in 4 dimensions, whose Penrose diagram is given in Fig. |], and is simply a 7r/2- 
rotation of the Penrose diagram for the Schwarzschild black hole. (The analysis 
for a charged brane can be done in analogy with the present case, with virtually 
no modification to the equations we present here, and with the same results.) The 
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Figure 1: Penrose diagram for the study of stability. 



metric for this spacetime is: 

ds 2 = - (l-^f) 1 dt 2 + (l-^f) dr 2 + t 2 d6 2 + t 2 sinh 2 ## 2 , (2.1) 
where 



t J V t 



h(t) = 1 - y • (2.2) 

The coordinate t is a time coordinate in regions I and III of Fig. [1], but is a spatial 
coordinate in the other two regions. 

The coordinates of eqs. (|2.1|) and (|2.2|) break down on the surfaces t = 2p, which 
correspond to the diagonal lines which form the boundaries of regions I and III. 
These are the horizons for this geometry. From the Penrose diagram it is easy to see 
that the line CDE is a Cauchy horizon, since the initial-value problem in region III 
does not uniquely determine the field evolution at points to the future of this line. 
Field evolution past the line CDE is not unique because it can also be influenced by 
signals from the time-like singularities. 

The Klein-Gordon equation for a massive scalar field propagating in the back- 
ground eq. (|2~T|) , is given by 

- ^-d M Ugg MN d N ] ^ + M 2 ^ = 

V 9 

in the time-dependent regions I and III. Now, since we are interested in the near 
horizon limit of the modes, it is convenient to write the KG equation in terms of 
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isotropic coordinates defined as r = t — 2p; the equation is then given by 



r 2 H + Vh 



d,. 



Vhh^d, 



^ + M 2 ^ = 0. (2.3) 



Here, for clarity, we denote hij(8, <fi) for the metric on the 2-dimensional maximally- 
symmetric hyperbolic space, and write gij(r,9,(p) = r 2 H + hij(9, 0) , where H + = 
1 + 2p/r. The relevant metric components are: g TT = —H + and g rr = H+ 1 . The 
functional form of the metric involved permits separation of variables, so we take 
ip(r, t, 9, 0) = e lcrr f{r) Lk(9, 0) , where a and k are separation constants determined 
by the eigenvalue equations: 



-die 



2 iar 



a 2 e lar and 



Vhh ij di 



k 2 L 



k- 



Both eigenvalue equations can be solved explicitly, and delta-function normalizability 
of the solutions require both a 2 > and k 2 > . The temporal eigenvalue equation 
then becomes: 



1 d 



r rrdf 



+ 



g"a 2 + 



A' 2 



T 2 H, 



+ M 2 



f = o. 



(2.4) 



Near the horizon, r — > and the asymptotic form is governed by the limits H + — > 

2p/r. The metric functions therefore reduce to g TT — > a T r~ 1 , g rr — > a r r and 
c<j — > The precise values of the constants ct r , a r and a u are not required, apart 
from the following ratio: ^ = r\. With these limits, the Klein-Gordon equation 
becomes, in the near- horizon limit: 



r 



2 1 



+ a T T 



M 2 + ^ 

a 2 , 



7 = 0, 



(2.5) 



If a 7^ 0, then the solutions are oscillatory, having the form / ~ r a ° , with ao = 
±2<Xa/ a T /a r . If a = 0, then a similar argument shows that the solutions are non- 
singular as r — > . 

One can now estimate whether an instability does exist by computing the energy, 
E = —u m d m i\) of the Klein-Gordon modes considered above, as seen by an observer 
whose velocity, u = Md t + Nd r , is well-behaved as it crosses the horizon. The 
normalization condition u 2 = — 1 in the vicinity of the horizon allows a determination 
of how M and N must behave as r — > (in isotropic coordinates) in order to remain 
non-singular. We find in this way u 2 ~ —a T M 2 r~ l + a T N 2 r , which is regular near 



r — > provided M 
then finds 



r 1 / 2 and iV 



-1/2 



near the horizon. With this choice, one 



—E = Md T ijj + Nd r *p ~ V ^ _1/2 



(2.6) 
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Using the asymptotic solution found below eq.( [2.5|) : ip ~ r a ° with a = ±«cta/ a T /a r , 
we see that E — > oo as the horizon is approached. This suggests that the stress- 
energy density of the mode under consideration diverges as well in this limit. As 
such, this mode is likely to destabilize the metric modes near the past horizon. 

We will now confirm this result by a full analysis of the metric perturbations. 
For later purposes it is convenient to define several quantities. We define the surface 
gravity of the solution to be 



1 

K = — 

2 

and the 'tortoise' coordinate as 



dh(t) 



dt 



(2-7) 

t=2p 



( - s /!r t+ i ln|t - 2p| - (2 - 8) 

where — oo < < +oo corresponds to the range 2p < t < oo. Notice that t* 
increases from past to future in region I, but decreases from past to future for region 
III. 

Our focus is on metric perturbations in region III of the Penrose diagram, which 
is to the past of the Cauchy horizon. We also focus on the 'axial' perturbations of 
the metric, which are defined as follows ||: Take one of the angular coordinates 
and change its line element by 

d<p 2 -> (d<p - q x dr - q 2 dt - q 3 d9) 2 (2.9) 

where the q^s are arbitrary functions of r, t and 9. The metric perturbations can 
be written, when specialized to the 'axial' modes, as a scalar equation for the field, 
<&(t,r,6) = t 2 h(t) (q 2> e — (fe^) sinh 3 9 (see || for details). Given the symmetries of 
the problem, the field equation for $ can be solved by separation of variables, with 
the field $(t,r, 9) decomposed as 

$(t,r,9) =r 1 Z(t)Q k (9)e i ° r , (2.10) 

with the functions Qk{9) and e tar , defined as solutions to the following eigenvalue 
equations II: 



d 2 ;„„ , , . , „ d ( 1 d t 

■>k 



jar = a 2 e iar ^ ^3 q q = j^q 



dr 2 ' ^ d9\smh s 9d9 

Writing the relevant equations in terms of the tortoise coordinate defined above, 
the temporal eigenvalue equation then becomes: 



dt 



2 
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where the potential in the previous equation is given by 

~2p 



2 h{t) + fc 



(2.12) 



Our interest is in the behavior of the mode $ near the horizon in region III 
(£* — > — oo), given its form in the asymptotic past (t* — > oo). We therefore need the 
asymptotic behavior of the potential in these limits, which is 

1/t 2 for i* — ► t — *■ oo 
e 2K °'* for — > -oo. 



T/^ \ 11 1UI ^ t ^ w , v 



Since the potential falls off faster than , the asymptotic behavior of the solutions 
of (|2.11| ) for t* — > ±oo is given by e ±lat * . 



Our initial condition must have only incoming waves in the asymptotic past, 
which for region III means for t — > +oo. Given this we wish to compute the coeffi- 
cients A(a) and B(a) which control the behavior of the solutions near the horizon, 
according to 

Z{U) -> e iat * U^t^oo (2.14) 

-> A(o) e~ icTU + B(a) e iaU U -»■ -oo . (2.15) 

(These boundary conditions look slightly odd compared to the Black Hole case, 
due to the reversal of roles between +oo and — oo which may be traced to the 
unconventional property that £ # becomes more negative into the future for region 
III (see Fig. 1). This difference of boundary conditions represents a key difference 
between the SO-brane geometry of interest here and other similar geometries, like the 
Reissner-Nordstrom metric.) 

In this way, the problem is reduced to that of scattering of waves on the potential 
of eq. ( 2.12j) . In order to distinguish transmission from reflection for a given initial 



wave we must distinguish the edges CD and DE of the Penrose diagram. In order to 
do so, we restore the r-dependence, e lcrr , of the initial wave, corresponding to a wave 
whose wavefronts initially move towards increasing r . It is convenient to express the 
solutions in terms of the null-like coordinates u = £* + r , v = £* — r , in which case 
the initial configuration is Z(t*,r) e mu for t —>■ oo, while the near-horizon limit 
of eq. fl2.15|) becomes 



Z{U, r) -> e iau + \B(a) - 1] e iau + A{a) e~ i(TV . (2.16) 

From the previous expression (or from looking at the figure) it is clear that the 
transmitted part of the wave will cross the edge CD , while the reflected part crosses 
edge DE. 

If we take a general, properly weighted, initial amplitude W(a) then the above 
formulae can be rewritten as 

Z(U,r)->X(v)+Y(u) (u,v->-oo) (2.17) 
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where 



X(v) = W(a)A(a)e- tav da, (2.18) 

J — oo 

/oo 
W(a) [B(a) - l]e iau da. (2.19) 
-oo 

To argue for instability we now compute the energy contained in the radiation as 
seen by a radially-moving inertial observer crossing the Cauchy horizon. The (n + 2) 
velocity, U, of such an observer is given by ||: 

u '=%=w Y u '--t-m lE2 - m1 ' 2 ' u '=°- (2 ' 20) 

where r is the proper time and we choose the negative sign for U 1 * in region (III) 
so that U is future-directed. Notice that the integration constant, E, which labels 
the observer's geodesic can be negative. 

A measure of the energy of the fluctuation $ as seen by this observer is given 
by T , defined by 

T = U"Z tll = U r Z :T . + U u Z, u , (2.21) 

or: 

l EZ * - - Kt)Y /2 Z,u] ■ (2-22) 
In terms of X and Y , in the near- horizon limit we have 

Z, r ^ X_ v + Y u , (2.23) 
Z, u -> ~X,- V + Y tU , (2.24) 



and so eq. ( |2.22|) becomes 



T - \ [X,- v (E + [E 2 - hf 2 ) +Y, U (E- [E 2 - hf 2 )] . (2.25) 

We now ask whether T diverges on one of the horizons, CD or DE. On CD, 
v remains finite, while u — > — oo. This means that, for E < 0, the term involving 
X jV of (p . 2 5|) remains finite while the term with Y )U could diverge. Hence: 

T CD -> -4p\E\Y tU e- K ° u (u^-oo on CD). (2.26) 

On DE, by contrast, it is u which remains finite, while v — > — oo. In this case, 
when E > 0, the term involving Y yU remains finite but the term with X „„ diverges. 
Hence 

Tbe -^4pEX_ v e~ K ° v (v^-oo on DE) . (2.27) 
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From these expressions we see that the existence of a divergence in T depends 
on the behavior of 

/oo 
W(a)iaA(a)e- mv da, (2.28) 
-oo 

/oo 
W(a)ia[B(a) -l]e ir7U da, (2.29) 
-oo 

near the horizons. The integrals may be evaluated by contour integration, with the 
contour closed in the lower half-plane for Y u (where our interest is in u — > — oo) and 
in the upper half plane for X_„ (where we care about v — > — oo). The result depends 
on the singularities of A(a) and B(a), which have analytic properties which can be 
quite generally determined using the arguments of Chandrasekhar and Hartle [§. || 
with virtually no modification. One finds in this way the function A(a) is analytic 
on the upper half plane except for the poles located at inKo , where n is a natural 
number, while the function B (a) is analytic on the entire upper half plane. This 
singularity structure is illustrated in Fig. 2. 



;k 



B(ct) A(o) 



Figure 2: Domains of analyticity of A(a) and B{a). 

Applying these results to the flux on surface DE, we evaluate the integral by 
deforming the contour into the upper half plane. The dominant contribution is then 
given by the closest pole to the real axis, which occurs at a = in . Then the integral 
giving X- v becomes 

oc e K0V . (2.30) 

We see, from ( |2.27| ), that !Fde, is therefore bounded as the observer crosses the 
horizon. 

Similarly evaluating the integral for Y^ u by closing the contour in the lower half- 
plane, we find contributions only from the real axis: ia = 0. Consequently, for W(a) 
analytic on the real axis, Y^ u is 0(1) near the horizon CD. It then follows from 
( |2.26|) that the flux Tqd at CD necessarily diverges. 

Although we have derived the instability for the lower Cauchy horizons, we also 
expect the same conclusion to hold for the upper horizons which bound region I of 
the Penrose diagram. We can argue this on grounds of continuity, as we follow the 
energy seen by a family of observers who cross the horizons in the vicinity of point 
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D of the Penrose diagram. We expect from this that the same infinite blue-shift 
seen by inertial observers for horizons CD and DE also extends to the other two 
horizons. 

This calculation confirms the preliminary work of ref. 0, which argued for infi- 
nite energy for Klein-Gordon modes on this metric. The present calculation shows 
that the same conclusion also holds for bona fide metric modes. The existence of 
this divergent energy strongly suggests that the horizon is unstable towards becom- 
ing a curvature singularity, due to the metric's back-reaction of these large energy 
densities: gravitational perturbations will introduce null-like singularities in the ge- 
ometry. This conclusion strongly resembles the same results for the instability of the 
RN black hole. 

The previous calculation can be extended to more general asymptotically-flat 
time-dependent backgrounds representing cosmological horizons. In general, the 
temporal eigenvalue equation that rules the perturbations will take the form of a 
Schrodinger equation like eq. ( |2.11|) . When the resulting potential is characterized 
by an asymptotic behavior given by (|2.13|) , the resulting cosmological horizon will be 
unstable under cosmological gravitational perturbations. Also, the analysis of 'polar' 
perturbations gives similar results as does the axial case presented above 0. 



3. Particle Production 



We now turn our attention to the time-dependent region (region I) which is to the 
future of all of the horizons. Because the metric in this region is not static, it should 
cause particle production for any quantum fields which propagate within it. We 
compute this particle production here, and show that it has features which resemble 
a thermal distribution whose temperature is given by the Hawking temperature, as 
defined in ref. || for the static part of the metric. We find in this way a connection 
between the properties of the time-dependent region I and the static regions which 
are separated from it by the horizons. These did not a priori need to be related, and 
such a relation seems even more odd if the surface dividing these regions represents 
a curvature singularity rather than a horizon. 

To this end we again consider the simplest non-trivial case, namely the SO-brane 
in four dimensions, with metric given by (|2.1|) , 

als 2 = - (l-y) 1 dt 2 + (l-y) dr 2 + t 2 d9 2 + t 2 sinh 2 ^0 2 . (3.1) 

Let us first recall the formal calculation of the Hawking temperature performed 
in ||. This requires the cancellation of a conical singularity, that can be obtained 
by requiring the proper periodicity for the Euclidean time in the Euclidean section 
of the metric's static regions. In the near-horizon limit the Euclidean metric for the 
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static regions takes the form: 



ds% = dR + KqR dr +••• (3.2) 

where R 2 = (2p — r) /2p and kq = l/(4p) denotes the surface gravity at the horizon, 
as in the previous section. Demanding no conical singularity at the horizon (R = 0) 
requires the Euclidean time coordinate r to be periodic r ~ r + 2n /k Qi leading to 
the Hawking temperature: 

Th = ^ = 5— ■ 3.3 

27T STCp 

Ref. || argued this temperature to be interpretable as the temperature of the particle 
distribution seen by static particle detectors in the static regions. 

We now compute a logically unrelated quantity: the particle production of a 
massless Klein-Gordon field caused by the time-dependent fields in region I. To this 
aim we consider the massless Klein-Gordon equation 1 : 

1 -A (V^gg^) $ = o, (3.4) 



which we again solve by separating variables to obtain the following mode functions: 
w mcr fc(t, r, 9, 0) = e l ( m< t> +(7r ) Q km {9) Vt{t) + c.c. . Here the integer m and real quantities 
a and k are the quantum numbers associated with the 0, 6 and r coordinates, 
respectively. We obtain the following 6 dependence: 

<d km (0) = a Qi iVm ^' 1] (coshtf) + b pi {VTTWI - l) (cosh#) , (3.5) 

where a and b are integration constants. P^(x) and Q r m {x) are the usual associated 
Legendre functions. In what follows we restrict our analysis to the simplest case, 
m = k = 0. 

The physics of interest lies in the t-dependent part, which can be usefully rewrit- 
ten as a Schrodinger-like equation by performing the substitution 



Q(t) = F(t) 



V 



y/t{t-2p) 



(3.6) 



It is also convenient to perform a change of independent variable x = t — 2p in order 
to place the horizon at x — 0. With these choices the t-dependent equation takes 
the form: 

d 2 F(r) 

- 1 ^-y^)F(x) = (3.7) 

with 

„, n. fa 2 (x + 2p) 2 2(x+p) 2(x+p)p 1 (x+p) 

- V ( x ) = Z3. + ZJZ-T^v?. + 



x 2 x(x + 2p) 2 x 2 (x + 2p) 2 x(x + 2p) x 2 (x + 2p) 2 

(3 



1 See irT[ for the analysis of the black hole case. 
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In order to solve this equation we replace V(x) with an approximate potential, 
V appTOX (x) , which is chosen to properly reproduce the asymptotic form of V(x) as x — > 
oo. (This is similar in spirit to replacing V(x) with one of its Pade approximants.) 
For these purposes we choose 



"^approx^) 



Ap Ap 2 , 2 

1 + — + -V I ° 

X x z 



(3.9) 



We compare the approximate potential with V(x) in Figs. 3 and 4, for different 
values of particle label a. As is clear from these figures, the approximate potential 
follows V(x) more closely the larger o is and the further x is chosen from the 
horizon [x = 0). Remarkably, even for op = 1 the potentials only deviate by 
a few percent right at the horizon, where the fractional deviation becomes (V — 
V appTaK )/V — > l/(16p 2 <7 2 + 1). For this reason we believe the approximate potential 
to more accurately capture the form of the Klein Gordon solutions near the horizon 
than would be possible using only an asymptotic expansion of the solutions in powers 
of \/{ap). 



Figure 3: Percent difference between the function and its approximation, op = 1 



Using the approximate potential the Klein Gordon equation becomes 
d 2 F(x) ( _ 4p 4p 2> 



dx 2 



1 + — + 

x x A 



o 2 F(x) 







(3.10) 



which can be solved exactly to give Whittaker functions [0 M x ^(z) and M. 



as the linearly-independent solutions. These are related to standard confluent hy- 
pergeometric functions according to 



(3.11) 
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Figure 4: Percent difference between the function and it's approximation, ap = 100 



The parameters x-, f 1 an d z are given in terms of p, a and x by the relations: 
X = —2iap, z = 2iax and /x = i/ij, with /i 7 = \ ^lQa 2 p 2 — 1 . Notice that x-, I 2 an d 
z are all pure imaginary so long as \ap\ > |, as we shall assume in what follows. 

For the purposes of a particle-production calculation we are interested in the 
combination of these functions which have positive and negative frequency near x = 
and x — > oo. It happens that it is M x _^{z) which is positive frequency near x = 
and M Xi(U (z) which is negative frequency, as may be seen from the small- z limit 



M X M = [1 + 0{z)\ . 



(3.12) 



The assignment of positive and negative frequencies follows once this singular 
part is re-expressed in terms of the time coordinate, x, in which case z^ + 2 oc 
xz exp [i/ij \og(x/p)] . (Recall the standard phase convention calls exp[—i<p(t)] positive- 
frequency when <f(t) increases with advancing time, such as for ip — cut with uj > 0.) 
For later purposes we also record here the useful identity: 



[M 



M_ 



e m z) 



(3.13) 



where we assume \ap\ > \ in order to use that all three of %, /i and z are pure 
imaginary. 

For x — > oo , on the other hand, it is the particular linear combination 



w x M = 



r(-2n) 



r(2/i) 



r(i/2- Ai - x ) — ' r(i/2 + / i- x ) 

which is positive frequency, as may be seen from its asymptotic form 

1" 



(3.14) 



W x M~z*e-'' 2 



l + O 



(3.15) 
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since e~ 2//2 = e~ l ^ x . That it is the absolute value of a which appears here follows 
from a careful treatment of the phase accumulated when z changes sign due to the 
branch cut at z = , using identities like eq. ( |3.13| ) . The additional phase associated 
with the factor z x does not change this conclusion. For instance, for op ^> 1 its 
effect is simply to change e~ l ^ x to e~ l ^ x * , where = x+2p \og{x/p) is the tortoise 
coordinate. 

With these preliminaries we may now proceed with the particle-production cal- 
culation. If we start with the mode expansion which is appropriate for large x, we 
have 

.(*,r,M)=E/**'['WW (*.r,M>H, (3 ' 16) 

m 

where dkma denotes the mode destruction operator and Ukmo(x, r, 9, 0) oc e mr W Xifl (z) . 
On the other hand, near the horizon we instead have 

«,,r^) = W**[ i ^,w(.,r,M) + c. ft 



(3.17) 



where bk ma are destruction operators and Vk ma (x,r,9,(j)) oc e mr M Xi _ M (z) . Particle 
production occurs because the expansion of Uk ma in terms of Vk ma implies that bk ma 
can be expressed as a linear combination of ak ma and a* k _ m _ CT . (It is conservation 
of m and a which permits only modes with opposite signs of m and a to mix in 
this way.) 

For simplicity it is convenient at this point to choose k = m = and to suppress 
the k and m labels. The decomposition of b a in terms of a a and a* a is found 
by choosing a particular o > and following those terms whose r-dependence is 
proportional to e tcrr . Keeping in mind that x{~ cr ) = — x( a ) = ~Xi z {~ a ) = 
—z(a) = —z and fi(— a) = /i(cr) = /i, and using eqs. ( |3 . 1 3| ) and ( ^.14| ) we find: 



xA~ z ) 
r(2/i) 



(3.18) 



r(| + M-x) 
r(2/i) 



a a M. 



z) + a! 



a --U-^J M_ Xi _ M (e- in z) 



z 2 



-Mr(l-2/*) 



(3.19) 



r(| + //-x) 

The last line gives the asymptotic form near x = . 

From these manipulations we see that the operators a a and b a are related to 
one another by: 



A(a) 



e 2 ^ a* 



(3.20) 



where A 2 = 1/[1 — exp(4i7r/i)] is determined from the normalization requirement 
[ba, b*] = [a a , a*] . Inverting this relation gives the expression 



A(a) 



b a + e 



(3.21) 
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which is the main result which is required for the particle-production calculation. 

We imagine preparing the field in the ground state as seen by observers crossing 
the horizon: b a \0) = 0, and then asking for the number of late-time particles which 
this state would contain. We find in this way our final result: 

(N a ) =(0|a> CT |0) = e4 J_ i , (3.22) 

where we recall /i 7 = | (16a 2 p 2 — l) 1 ^ 2 . 

This result is very suggestive of a thermal form. Indeed if we write uj = [a 2 — 
l/(4p) 2 ] 1//2 , so uj ~ o when up ^> 1, then eq. ( |3.22| ) is precisely thermal, 

if uj is interpreted as the particle energy. (This interpretation is natural near the 
horizon where v a ~ ~ e -%niiogx ^ e -2ipujiogx ^ e -iwx* ghowg that uj is the 

eigenvalue of the operator idt„.) Since our derivation assumes /i is pure imaginary 
it breaks down for \ap\ < 4, where uj becomes imaginary. Our approximate form, 
Kpprox, also provides a worse description of the full result, V, for a this small. 
Consequently we cannot yet say whether these modes are also pair produced at late 
times. 

Even more remarkably, the corresponding temperature is 

T = -L , (3.24) 
8np 

which is exactly the same result obtained earlier by euclideanizing the metric in the 
static regions. A priori these did not have to agree since the euclidean calculation 
describes the particles seen by an accelerating, static observer behind the horizons, 
while in the present instance the temperature corresponds to the distribution of 
particles which are produced by the time-dependent fields in region I. 

4. Discussion 

The metrics studied here were proposed in ref. with an eye to using their time 
dependence for cosmological applications. Clearly, the classical instability of the 
horizon towards singularity formation diminishes the cosmological impact of these 
solutions. In particular it prevents the passage from region III (past, contracting 
time-dependent solution) into region I (future, expanding time-dependent solution) 
in a way which does not hit a singularity. This result is consistent with the strong 
cosmic censorship conjecture, since the observer crossing the Cauchy horizon, decou- 
pling from her past history, would otherwise find a naked singularity, from which 
information can come, and yet would be able to avoid the singularity. It would be 



interesting to see if there are cases that avoid this problem |T3, WA 
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A similarly wet blanket is thrown on the S-brane interpretation of this geometry 
in terms of a rolling tachyon field [[l], [15] . In this interpretation we imagine the rolling 
of the tachyon field from one minimum of the potential at t — > — oo to the other 
minimum at t — > oo . Each vacuum could then be identified with the asymptotic, flat 
infinite past and future of Fig. 1 respectively. The local maximum of the potential 
would then be identified with the horizon. Our result presents an obstruction to this 
realization of the rolling by making it impossible to miss a singularity in between. 
(This singularity problem is also shared by other proposed S-brane geometries.) 

The existence of the singularity need not invalidate the interpretation of region 
I as describing the metric produced by the late-time rolling of a tachyon from a local 
maximum to a later local minimum, however, since this part can be described purely 
by the geometry in the future of these singularities. Following this interpretation, our 
result for the particle production could be relevant to the determination of particle 
production after tachyon condensation. This would be particularly interesting for the 
string scenarios of hybrid inflation from D-brane interactions, as proposed in pi . In 



this case the particle production could lead to the determination of the re-heating 
after inflation (for a recent discussion of reheating from tachyon condensation see 

Q)- 

We find our particle-production result to be intriguing in its own right, due 
to the thermal character of the produced particles, and the connection which it 
indicates between the temperature of this distribution and the temperature obtained 
by euclideanizing the metric's static region. This connection is all the more intriguing 
given the classical instability which we find, which is likely to convert the intervening 
horizons into curvature singularities. 
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